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"v£>" \ Abstract 

I investigate the evolution of finite temperature, classical Yang-Mills field equations un- 
der the influence of a chemical potential for Chern Simons number Ncs- The rate of Ncs 
diffusion, T d , and the linear response of Ncs to a chemical potential, T M , are both computed; 
1—5 ■ the relation = 2T^ is satisfied numerically and the results agree with the recent mea- 
t— i . surement of by Ambjorn and Krasnitz. The response of Ncs under chemical potential 
^ \ remains linear at least to fi — 6T, which is impossible if there is a free energy barrier to the 
| motion of Ncs- The possibility that the result depends on lattice artefacts via hard thermal 
loops is investigated by changing the lattice action and by examining elongated rectangular 
lattices; provided that the lattice is fine enough, the result is weakly if at all dependent on the 
specifics of the cutoff. I also compare SU(2) with SU(3) and find Psu(3) ~ 7(a s /a w ) A Tsv(2)- 
§ ! PACS numbers: 11.10. Wx, 11.15.Ha, 11.15.Kc, 11.30.Fs 
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p; 1 Introduction 

cu : 

Baryon number is violated in the Standard Model 0. While it is conserved to all orders 
. ^ | in perturbation theory, nonperturbative effects involving topologically nontrivial gauge and 
Higgs field configurations permit its violation at a rate which, at zero temperature, is sup- 
pressed by an exponent of order exp(— 4n/aw)- While this rate is far too low to have any 
phenomenologically interesting consequences, it should be much larger at high temperatures, 
perhaps high enough to explain the matter abundance in the universe. 

I will briefly review baryon number violation in the Standard Model. Because of the axial 
anomaly, the Baryon number current is not conserved, but satisfies 

d^B = Z §^e^F^, (1) 

where Np is the number of generations of fermions and F is the SU(2) field strength tensor. I 
have left out the hypercharge fields, which are irrelevant because they do not have nontrivial 
topological properties. For smooth field configurations the righthand side of this equation is 
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a total divergence, 
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The charge associated with this quantity, / d?x(Ko), is called the Chern-Simons number, 
A^cs- It possesses the important property that, if the field configuration is pure gauge at 
two times t — to and t = t\, then 

N CS (h) - Ncsfo) = £ dt J d'x-^e^F^ (3) 

is an integer (assuming only that the gauge field is everywhere smooth). Ncs is invariant 
under small gauge transformations (those which can be built up infinitesimally) and changes 
by an integer under large gauge transformations; the difference between Ncs a ^ two times is 
gauge invariant, because the righthand side of the last equation is. 

The action of a Euclidean process which changes Ncs is at least / (1/4)F" i 7 ^ > 8ii 2 /g 2 , 
and is in fact even larger, because when is nonzero, the Higgs field gradient generally 
isn't either. If the configuration is of greater spatial extent than the inverse W mass then the 
Higgs field gradients dominate the action, which will be 3> 8n 2 /g 2 . For configurations of very 
small spatial extent, the asymptotic freedom of the weak coupling g again makes the action 
grow. Therefore, unlike QCD, topology change in the electroweak theory is exponentially 
suppressed, as mentioned above. 

However, at finite temperature it is no longer relevant how large the action of a configu- 
ration is; the system is excited and there is energy available to make topological transitions. 
In the broken electroweak phase there is a free energy barrier to such transitions 0, [| , which 
proceed at an exponentially suppressed rate determined by the temperature and Higgs mass. 
At higher temperatures electroweak symmetry is restored, so the Higgs gradients no longer 
inhibit topology changing transitions. For sufficiently large configurations, with spatial ex- 
tent 0(l/g 2 T), the free energy barrier, if any, is parametrically order unity, and baryon 
number may be violated readily [[| . This is relevant because any B + L (Baryon plus Lepton 
number) present in the early universe would then be erased, and, during the electroweak 
phase transition, any B + L separation by the wall of a broken phase bubble would produce 
a baryon number asymmetry, because the excess on the symmetric phase side would be 
destroyed, while that on the broken phase side would be preserved^], [|. 

The rate which is relevant to this baryogenesis scenario is probably the rate at which a 
chemical potential for baryon number, caused by the separation of baryon number across 
the wall, induces baryon number changing transitions^. That is, we need to know 



2 only probably, because in the realistic problem the chemical potential may extend over a region of space 
or time small enough that the large volume or large time limit are not attained 
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where by () M I mean the expectation value when there is a chemical potential \x for Ncs 
(equal to Np times the chemical potential for baryon number, by Eq. (|T])), and V is the 
volume over which the chemical potential exists. 

A detailed balance argument shows that this rate is half the rate for Ncs diffusion 0, 

rj . lm ««h«f) (5) 

t^oc Vt 

in the absence of a chemical potential. The argument is quite simple in the broken phase; 
the general case is reviewed in section 2. The nice thing about Td is that it is straightforward 
to measure it for the classical bosonic theory at finite temperature by lattice techniques, and 
there is reason to believe that the behavior of this classical system corresponds well to the 
infrared dynamics of the actual finite temperature, quantum Yang-Mills Higgs system||. 
Based on these ideas, and following previous work||, Ambjorn and Krasnitz have re- 



cently performed extensive numerical simulations of the lattice Yang-Mills system |L0| and 
have demonstrated convincingly that there is a lattice spacing independent classical value 
for T d . On parametric grounds, T d = K d (a w T) 4 ; they find n d = 1-09 ± 0.04. This work 
apparently answers the question relevant to baryogenesis, but it leaves some open questions. 
For instance, because the theory has linear divergences, the lattice spacing independence 
does not necessarily assure the absence of lattice artefacts. Also, a physical explanation of 
what sets the rate is lacking. It is smaller than the analytic estimate found in /citeFarrar- 
Shap by a factor of 40, and in comparison to the length 21 of a cubic lattice large enough to 
eliminate finite volume corrections in its computation, the rate is about Z~ 4 /100 (I include 



the 2 to account for periodic boundary conditions) It is possible but not certain that 



this apparent suppression arises from a (parametrically order unity) free energy barrier to 
winding number change, as proposed for instance in |TjJ . 

An alternative technique, which can be used to shed some light on these matters as well 
as to check the result of Ambjorn and Krasnitz, is to apply a chemical potential for Ncs 
on the lattice and look at the rate at which Ncs drifts. The method encounters some fairly 
serious technical difficulties M, which I will discuss at some length in section 3; I present a 
resolution, which the numerical results of section 4 demonstrate to be successful. I will also 
show that there is in fact no free energy barrier to topology change in Yang-Mills theory at 
finite temperature. Section 5 will be devoted to investigations of possible ultraviolet lattice 
artefacts in the classical lattice simultations. In section 6 I extend the technique to find the 
strong N C s diffusion rate (I replace the group SU(2) with SU(3)). In section 7 I conclude. 

1 present a simple, efficient thermalization algorithm which works for both the SU(2) and 
SU(3) cases in Appendix A. 

2 The Two Rates 

As mentioned above, there is a general relation between the Ncs diffusion rate T d and the 
linear response rate to a chemical potential, T^. This relation was demonstrated already in 
0, and although their argument explicitly assumes that there is a substantial free energy 
barrier, or, in the symmetric phase, that the rate is very small, that assumption is actually 



immaterial to their argument, as has been shown more recently in where the relation is 
derived in full generality. 
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Here I will present an alternate proof, which displays that the relation between the two 
rates only relies on the system's periodicity in Ncs and its ergodicity. 

Let us assume that the generalized coordinate Ncs is i n contact with a large thermal reser- 
voir, which evolves ergodically, so that for large t, {(Ncs(2t) — Ncs (t)) (Ncs (t) — N C s(0))) = 
0(1) does not grow as a positive power of t. This is the condition that the motion of N C s 
does not have arbitrarily large time correlation and is necessary for a diffusion rate to be well 
defined; in fact it immediately follows that ((Ncs(t) — Ncs(Q)) 2 ) = 0(t). I also assume that, 
for any operator O, (Ncs(t)O(0)) = 0(1) does not grow as a positive power of t. (These 
conditions should generically be satisfied for an ergodic system.) In this case, the probability 
for a system to go from N C s(0) at time to N C s(t) at time t depends at leading order for 
large t only on those values and t, and the probability that the Chern-Simons number is 
Ncs(t) a t time t is, in terms of the probability distribution V(N C s(0), 0) for its value at 
time 0, 

V(N cs (t),t) = [ dNcs(0)V(N C s(0),0)G t (Ncs(0),Ncs(t)) , (6) 



where Gt(Ncs(Q), Ncs(t)) is the probability that the system, initially at iVcs(O), will arrive 
at time t at Ncs(t)- Gt is invariant to a simultaneous shift of both arguments by an integer, 
and in the absence of a chemical potential, it is invariant on exchange of its arguments (by 
time reversal invariance), and, when one argument is 0, is invariant under a sign change in 
the other (by parity invariance). 

To determine the rate of response to a chemical potential, I want to consider the equilib- 
rium probability distribution when a term fiNcs is added to the Hamiltonian. Unfortunately 
the energy is then unbounded from below and no equilibrium probability distribution exists. 
I will therefore add a term eN cs to the Hamiltonian so that the theory will be well defined; 
in fact, in the real theory in finite volume there is such a term, because, as Ncs changes, 
the number of fermions slowly builds up and changes the chemical potential. I take e to be 
much less than any other quantity in the problem and consider the equilibrium probability 
distribution V^Ncs)- The chemical potential changes the Boltzmann factors for different 
states, so in relation to the distribution without chemical potential, 

^( N csi) Vo(Ncsi) f (N N Wp -/^/v N \ 

r> l at \ - r> l AT T J periodic [NcSl, N C S2> X exp — {1\ C S1 ~ N C S2) 

rn{^CS2) ro{^CS2) J- 

Vo(Ncsi) (1 _ fJ.(N CS i - N CS 2) V 2 (N CS i - N CS2 ) 2 

^ + ^fi e no A i C (Ncsi,Ncs2)^ , (7) 

where the function / pcr i dic is periodic on integer changes of either argument and equals 1 
when the arguments coincide. Its expansion only contains terms even in /i. 
Also, the hopping probability G t takes on \i dependent corrections 

Gt* = GU ±G\ + ^- 2 G\ (8) 



where G\ satisfies 



Gl(0,N C s(t)) = -Gl(0,-N cs (t)) (9) 
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because the chemical potential is odd in Ncs- There is also an order e correction which will 
not be important here. 

Using the condition that, for the equilibrium distribution (with or without //), V(Ncs) 
is time independent, I find 

^(0,0) = 7yo,t) = I dN C s(0)V^Ncs(0),0)G t ^N cs (0),0), 

V (0,0) = V {0,t) = JdN C s(0)Vo(N cs (0),0)G t (N C s(0),0). (10) 

Dividing through by V(0, 0), I get the equality 

/ ^'Sy' ^* ° )= I dN cdO) M ^ 0) G°(Ncs(0), 0) , (11) 

which I can expand in powers of /x, using Eqs. (|7]) and (|J). 
At order /i, 

I^o^^^Nos^J^os^kfi^O), (12) 

which follows automatically because both sides are zero; the first because, in the theory 
without chemical potential, {Ncs) — 0, and the second because G\ is asymmetric. 
At the next order, 

S iN - v -w N -f G ° M) = /* ! w w!fcl) (13) 

- / dNcsf^ J^^ G'iNcs, 0) (14) 

The term on the lefthand side is half the mean squared change in Ncs between time and 
time t, in the system without chemical potential, which is VtTd/2; the righthand side is 
VtT/j,. The first line is the desired relation, 1^/2 = T^. All that must be shown is that the 
remaining terms do not grow linearly with large t. 

For (|TJD this holds because fp eviodic V (N C s, 0)/V Q (0, 0) is bounded, and / dN C sG t (N C s, 0) 
is 1 because it is the probability that if the system starts at = 0, it will end with any 
value of Ncs- Hence the term as a whole is bounded, independent of t, and cannot grow 
linearly in the large t limit. 

For the second term, (|T5D, we can see that it is bounded at large time by considering the 
evolution of the system with chemical potential, but when the initial probability distribution 
is Vq, the equilibrium distribution for zero chemical potential. Since the distribution starts 
out periodic, it must remain periodic, and 

P(0,t)/P(0,0) = J dNcs ^^ (G° t (Ncs,0) + ±G\{N oa J) + ^G 2 t (N cs ,0)) (16) 
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must remain bounded. The G\ term drops because the initial distribution is symmetric, and 
the term gives 1; the G 2 t term is the one we want to know, and we see that it is bounded. 
In fact, I expect that the probability distribution will approach Vo(Ncs) * /periodic at large t, 
so the two extra terms (|HD and ( p~5D probably cancel; but it is sufficient to show that they 
remain bounded at large t, which I have now done. Hence the relation Td = 2T fl follows. 

Note that it is essential in the above argument that the typical distance that Ncs moves, 
over a time long compared to the thermalization time, is much less than T/fi, a condition 
on the size of \i. If there is a free energy barrier between integer values of Ncs, then the 
important typical distance which Ncs moves will be at least 1, and we can only be assured 
of linear response to /i for // < T. 

With a few additional assumptions we can explore the behavior of the system for larger 
/i, when there is a large free energy barrier. I will assume that the free energy barrier is 
significantly larger than T, and that the motion of Ncs is quite heavily damped, so that if 
one barrier is breached, enough of the energy used (and, at nonzero /i, gathered) in breaching 
it is lost that the system is unlikely to breach the next barrier before becoming thermalized 
in the new local minimum. In this case, the probability density will become equilibrated 
about the local minimum, and the probability to go over a barrier will be a kinetic prefactor, 
times the Boltzmann suppression to get to the top of the barrier. If /x is a small fraction of 
the barrier height, then the shape of the barrier will not be modified much, and the kinetic 
prefactor will be approximately equal to its equilibrium value. In this case, the rate of Ncs 
change will be 

(Ncs) = V(T imd - r backward ) ~ -L- (e# - e&) = -T d V sinh A (17) 

which of course has the correct small \x behavior. However, as /i becomes on order or greater 
than T, the rate rises, eventually greatly exceeding the linear extrapolation. The behavior 
will depart from the sinh function when the energy gained from the chemical potential in 
hopping a barrier becomes comparable to the loss from friction, or when the shape of the 
barrier is significantly modified by the chemical potential term. Note that if, contrary to my 
assumption, the friction is weak, then most barrier jumping events are followed by another, 
and the diffusive motion will have very long time correlations. This is not observed in 
the simulations. I will just comment that, in this situation, the very large /i behavior of 
the system gives a rate which differs from the linear extrapolation by on order exp (barrier 
height). 

What happens when there is no barrier? I will mention the simplest possible model; that 
Ncs moves in a flat potential under a linear frictional force and a random force to keep it 
thermalized, 

N cs = -lN cs + f (18) 

where / is independent of N C s and has (f(t)f(0)) = for t ^ 0. Its magnitude is chosen 
so that Ncs w iH De thermalized at the correct temperature. Adding a chemical potential 
produces a term —[iNcs on the righthand side. Averaging over realizations of /, the terms 
Ncs and / drop out of the determination of (Ncs), which is strictly linear in /j. This 
behavior will eventually break down only because a realistic / does have a correlation time 
and a realistic friction is somewhat nonlinear and has some memory. 
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Now I will briefly discuss evaluating T^. It is a linear response coefficient, and as such 
we might expect that it can be computed from a correlator of the equilibrium theory. In 
fact this is the case (for the continuum theory); by considering modifying the Hamiltonian 
by H — > H + fiNcs and applying the techniques of Kubo |16[], we find 

r, = lim T ({Ncs(t),Ncs(0)}) 

where {, } are the Poisson brackets, 

M m _ y 9A dB dB_JA_ 

1 ' 1 Y% (0)^(0) d qi (0)d Pi (0)' 

The derivatives in Eq. (|T9|) are with respect to the coordinates and momenta at zero time. 
The poisson bracket we need then includes expressions such as dNcs{t)/dp(0). To evaluate 
numerically the contribution to this expression from one pair of conjugate variables in one 
representative configuration, we would have to make a small change to one momentum and 
then evolve the system forward by time t, and compare Ncs to its value if we had performed 
the evolution without changing the momentum. Since we must sum over all coordinates and 
momenta (and on a lattice there are a great many) and over an ensamble of initial conditions, 
this will not be a very efficient way of measuring T^. A more practical way would be to 
perform all the terms in the sum on coordinates at once, by adjusting each coordinate and 
momentum according to the derivative of Ncs with respect to the corresponding momentum 
or coordinate. Provided that we make small changes, the final Ncs should depend linearly 
on each change, contributions arising because of the product of two of the changes we made 
being suppressed by a power of how small the changes were. We can also rely on the ergodicity 
of the system and apply the adjustments to coordinates and momenta at each time step, 
rather than at only one. This appears to be the only efficient way of evaluating Eq. (|19|); 
it is, however, equivalent to evaluating the evolution of the system under the action of the 
chemical potential. So we gain nothing by considering equilibrium correlators, and the best 
way to proceed numerically is to evolve the system with a chemical potential term added to 
the Hamiltonian. 



3 The Model 

Now I turn to the problem of how to measure T M . I will only be concerned with finding the 
high temperature limit. Far above the electroweak phase transition temperature the Higgs 
boson takes on a substantial plasma mass and probably does not influence the evolution of 
very infrared modes, so it will be sufficient to evolve the classical Yang-Mills field equations 
under the influence of a chemical potential for Ncs- 

To prepare for the task I review the method of evolving the Yang-Mills equations when 
there is no chemical potential. The evolution of the fields is not completely specified by 
the field equations because they respect the freedom to choose a gauge; specifying a gauge 
renders the evolution unique. The most convenient choice is the gauge which will make A 
everywhere zero; in terms of the gauge fields A 1 and their conjugate momenta, the electric 
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fields E\ and suppressing group indicies, the continuum equations of motion in this gauge 
are 



~dt 

dE l 

~~dT 



E' 

-D J F ji (21) 



where F is the field strength tensor, and I have written the space components in terms of 
a positive three dimensional metric. In addition, the time component of the Yang-Mills 
equations of motion, which arises from minimizing with respect to variations in A , enforces 
a first-class constraint, 

D i E i = . (22) 

On the constraint surface, this Gauss constraint commutes with the equations of motion, so 
it need only be enforced on the initial conditions. However, it is important that this property 
be preserved numerically when we actually evolve the field equations. 

It is also well known how to make a lattice version of the Yang-Mills field equations 
without a chemical potential fl7 |. The gauge fields become the link matricies U l , and the 



electric fields are adjoint fields which I will take to lie at the beginning of each link, so the 
link update rule is 

tP = B^iTjP . (23) 

(The left index of U lies at the beginning or basepoint of the link and the right index lies at 
the end. To parallel transport Ei to the end of its link, we then take UjEiUi.) 
The Kogut-Susskind Hamiltonian is 

tf = E^ + El-^Trf/ D (24) 

where the first sum is over all links and the second is over all elementary plaquettes, and 
U a means the product of the link matricies running around the plaquette. The coupling 
constant is absorbed into the inverse lattice temperature (3l which will be used to thermalize 
the system; with these conventions (3l corresponds to the continuum value (3l = 4/(g 2 aT), 
a the lattice spacing. 

The Gauss constraint at each lattice point is 

£ Ek{+) -J5fc(-) = (25) 

i 

where Ei(+) is the electric field on the link running forward out of the site and Ei(-) is 
on the link running into the site from behind, parallel transported to the lattice site. This 
linear combination of E fields generates a gauge transformation of the U fields at the lattice 
point, and the requirement that it vanish can be understood as arising from our using up our 
freedom to make time dependent gauge transformations in choosing the temporal gauge. A 
simple leapfrog algorith for the Kogut-Susskind Hamiltonian, as used in |9|, |IU| , identically 
preserves the Gauss constraint. Krasnitz has recently developed an algorithm, based on 
a set of Langevin equations, for thermalizing this system to inverse temperature (3l while 



identically preserving the Gauss constraint; for details see [JTT . 
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To extract the value for Td from the realtime evolution, one must keep track of the change 
in Ncs, which has a gauge invariant continuum definition, given in Eq. (^). In the leapfrog 
algorithm, where the U{ are defined at integer multiples of a time stepsize At and the electric 
fields live at half integer values, the simplest lattice realization which is invariant under the 
lattice point group is 

AJV CS _ 1 / Bf(t + At/2) + Ef(t - At/2) 11 \ 

~~^k { 2 5 k 2 " °) ' (26) 

where 8Dj means a sum over the 8 plaquettes which begin and end at an endpoint of the 
link and run in the plane perpendicular to the link, see figure |I|. (The factors of 2 differ from 
the continuum version because the lattice quantities are interms of r, while continuum fields 
are defined in terms of r/2.) 

Now I will discuss the relationship of this model to the physical, quantum system we hope 
it will simulate. The degrees of freedom of the thermal, quantum Yang-Mills system fall into 
3 categories. There are thermal energy particles with momenta characterized by the scale T, 
which interact very weakly with the other excitations, carry almost all of the energy of the 
system, and have ultrarelativistic dispersion relations. There are more infrared fields with 
momenta on order the plasma frequency ~ gT, which have substantial occupation numbers 
but interact perturbatively with the other degrees of freedom, except that forward scattering 
with the thermal particles make order unity modifications to their dispersion relations. And 
finally there are very infrared modes, with momentum characterized by the scale g 2 T, with 
large occupation numbers and fully nonperturbative mutual interactions. 

On the lattice there are also three characteristic scales: excitations with wavelength order 
the lattice spacing are weakly coupled, contain almost all of the system energy, and travel 
under the lattice dispersion relations; excitations of wavelength order a/fli, interact with 
each other nonperturbatively; and the intermediate scale again has perturbative interactions 
but dispersion relations substantially modified by interactions with the shortest wavelength 
modes. The motion of Ncs i n the symmetric electroweak phase depends on the most infrared 
modes. The idea of Grigoriev and Rubakov is that in the quantum system these have such 
high occupation numbers that they should behave approximately as classical fields, and 
should therefore be correctly modeled by the classical lattice system ||. If the high energy 
modes are only important as a thermal bath and for their corrections to dispersion relations, 
then the substitution of classical lattice modes for quantum continuum ones should not 
matter much. 

Ambjorn and Krasnitz have also noted that the thermodynamics of the lattice model 
bears close resemblance to that of the full quantum system. In particular, if one introduces 
Lagrange multipliers A to enforce the Gauss constraints and then performs the (Gaussian) 
integration over the electric fields, the partition function is almost the same as the partition 



function of the full theory, in the approximation of dimensional reduction [|T0| , |18f . The 
difference is that the classical, lattice theory has zero bare Debye mass. However, a Debye 
mass, linearly divergent in the lattice cutoff, is induced by the high energy excitations. Its 
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value is 

m 2 D (lattice) = Ag A T S = 3.1759114. (27) 

A-KCL 

The lattice system therefore does in fact reproduce the thermodynamics of the full theory, 
for the right value of the lattice constant aQ. Yang-Mills Higgs theory has m D = 5g 2 T 2 /6, 
obtained in classical lattice Yang-Mills theory by (3l — 7.8. The full standard model has 
m 2 D = llg 2 T 2 /6, corresponding to (3l — 17.2. The independence of the rate on (3l 
therefore indicates that Yd does not depend on the Debye mass. I will return to this point in 
section 5. I will also discuss a potential problem, involving the functional form of the "hard 
thermal loops" induced by the high energy modes. 

Now I turn to the problem of adding a chemical potential for Ncs- For the continuum 
equations of motion this is straightforward; one adds a term fiNcs to the Hamiltonian, and 
the evolution of A, E are modified by 

0)+ f i{A i ,N cs } = E i + 0fi 

0) + fi{E\ N cs } = D J I '" - -^-e ijk F jk = -D j F ji - -^-B i (28) 

where B l = (l/2)e ijk F jk is the magnetic field in the i direction. Fortunately, for smooth 
fields 

^ k D*F jk = (29) 

by the Bianchi identity, and so the additional term still preserves the Gauss constraint. 

An equivalent approach is to add to the action a term 9{t)e^ il/a pF^ F" 13 . If 6 were a 
constant, then this term would be a total derivative and would not change the equations 
of motion; but if it varies in time, then when we integrate by parts in time while deriving 
equations of motion the time derivative acts on 9, generating the B l term in the equations 
of motion. Proceeding in this fashion also makes it clear why the added term still commutes 
with the Gauss constraint; we should rederive the Gauss constraint with this added term in 
the action, but because the index of the antisymmetric tensor is used up by A , no time 
derivatives arise to act on 9, so the Gauss constraint is unchanged. 

When we try to implement a chemical potential for Ncs on the lattice, we encounter 
trouble. The most obvious way is to find a lattice definition of N C s and add it to the 
Hamiltonian. However, there is no lattice definition of Ncs which gives the desired behavior 
under gauge transformations. Ncs should be invariant under small gauge transformations 
and change by an integer under large ones; but on a lattice, a gauge transformation is a choice 
of one member of the gauge group at each lattice site, and since there are finitely many lattice 
sites and the gauge group SU(2) is path connected, any gauge transformation can be built 
up infinitesimally. Therefore, either Ncs must be a constant, or it must vary continuously 
under gauge transformations. Also, we would like Ncs to change by the same amount 

3 In that paper the coefficient is 5 rather than 4, because they work in Yang-Mills Higgs theory, and the 
interaction between Aq and Higgs fields contributes 1 to the coefficient. 

4 This statement is not quite correct, because the lattice theory also contains power law divergences in the 
energy density and in the expectation values of composite operators; this point is discussed in |is|| , which 
also shows that these divergences do not interfere with the extraction of physically interesting quantities. 



dA l 

~dT 

dE i 

~df 



dA l 

~dt 

dE 1 

~dt 
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on each application of the same gauge transformation; but any lattice Ncs is a function 
of finitely many variables, the Ui, which live on compact manifolds, and any continuous 
function on a compact domain is bounded. If it changed by some constant amount under a 
gauge transformation, then repeated application of the gauge transformation would drive it 
to infinity, a contradiction. 

A second idea for implementing the chemical potential on the lattice pOf is to use this 
definition of the lattice e^^F^F^ and the alternative implementation of the chemical 
potential in which we add 6(t)e >iual3 F^ v F aj3 to the (lattice) action. The Gauss law, which is 
derived from the lattice action, is modified in a way such that it will be preserved by the 
new equations of motion. 

This technique suffers from a new problem on the lattice, which is that the lattice def- 
inition of e^^F^Fap, Eq. (p6|) , is not a total derivative. The simplest counterexample 
configuration contains just three non-identity link matricies; two, with values 1 + eiri and 
1 + eir 2 , lie on the two most removed links of a plaquette perperndicular to some link, and 
the third, with value 1 + eir 3 , lies along the link, backwards one time step (see figure ||). 
The space integral of e^ 1101 ^ F^ y F a p for this configuration is nonzero even though all links at 
large distances are the unit link. The definition is therefore not a total derivative^. This 
is a serious problem, because it means that even if 9 is constant, the Oe^^F^F^ term 
will change the dynamics of the system; when 9 varies with time, then at late times it will 
become large, and these (spurious) changes will actually dominate the dynamics. 

It is not surprising that this technique did not work. The argument in the continuum 
theory that e^ uaP F^F a p is a total derivative relies on the fields being smooth, a concept 
which is lost on the lattice, and its failure to be a total derivative explains why Ncs is not 
well defined. Recall that the chemical potential for Ncs emerged anyway by integrating out 
chiral fermions at nonzero chemical potential, and that it is impossible to implement chiral 
fermions on a lattice, a fact which may be related to our problems here. 

Another idea, considered in 0, is to abandon the hope of deriving equations of motion 
from a Hamiltonian and to try to find a lattice implementation of the continuum equations 
of motion, Eq. P%|). A natural choice for Bi is the average of the eight plaquettes used in 
Eq. (p6|) ; the equations of motion are then 

U%t + At) = exp(iAtr a ^)C/ l '(t) 
K(t + At/2) = El{t - At/2) - £ (\^{-ir a )U a ) - J±- £ (\^{-ir a )U n ) .(30) 

4n . \Z / 107T gn . \Z / 

Here X^n^ means a sum over the 4 plaquettes which contain the link i, with orientation so 
as to contain U and not l/J , and SsDj nas the same meaning as previously. This technique 
is nice in that it corresponds to the physical meaning of a chemical potential term, that the 
E fields should be modified in accordance with the B fields so that the energy in the fields, 
E ■ E/2, is modified by E ■ 5E = -/lAtE ■ B/(2ir 2 ) = -fxAN CS - 

5 For any definition of the space integral of e^ ual3 F^Fup which consists of a sum over lattice sites or links 
of some local gauge invariant operator, a similar counterexample can always be found; find three links which 
all appear in the evaluation of one, but only one, site or link, and shift their link matricies from the identity 
using three orthogonal Lie algebra elements; then they will contribute to J FF at this point but not at any 
other, and the definition will not be a total derivative. Note that the nonabelian nature of SU(2) is essential. 
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The problem with this plan is that the evolution does not preserve the Gauss constraints, 
so it excites unphysical modes. Again, the nonabelian nature of the theory is essential; in 
the abelian theory, the lattice definition of e^ va ^ F^Fap is a total divergence, and the space 
divergence of the magnetic field at some point, which is the sum of all plaquettes on the 
surface of a box one lattice spacing around the point (see figure |]), is the boundary of a 
boundary, and vanishes identically, preserving the Gauss constraint. But in the nonabelian 
theory, while -Dj-B; is still the sum of all plaquettes around a cube, there are commutator 
terms which spoil the cancellation. That is, B l a is contaminated with nonrenormalizeable 
operators such as a 2 D l D l B l , and D l D l D l B l need not vanish, since [D l ,D l ] does not. Since 
the extra terms are nonrenormalizeable, they will have vanishing influence on the infrared 
physics as the lattice spacing is made smaller; but if they introduce ultraviolet divergences 
in the measured value of then they will spoil the calculation. 

If we ignore the violation of the Gauss constraint, or if we (dissipatively) remove the 
Gauss constraint violation occasionally or incompletely, we risk measuring changes in Ncs 
arising from the excitation of unphysical modes as well as from genuine topology change. 
The problem is that the combination of electric fields which produce a Gauss constraint 
violation has essetially no restoring force. If the chemical potential term causes this E field 
to change, then E ■ B will be nonzero for the linear combination of E which make up the 
Gauss constraint; if there were a restoring force then E would oscillate and the time integral 
of E ■ B would vanish, but instead E will rise and the integral need not vanish. If the Gauss 
constraint is dissipatively enforced occasionally, it is quite easy numerically to establish 
how much of any accumulated Ncs corresponded to this process. Because of the chemical 
potential, the system energy shifts by liAN C s- Some of this energy goes into exciting Gauss 
constraint violations; the amount is /i times the amount of ANcs which arose from exciting 
the unphysical modes. If the Gauss constraint violations are quenched, then the loss of 
energy in the quenching corresponds to the amount of spurious ANcs- Numerically I find 
that this amount is a substantial share, if the Gauss constraint violation is only quenched 
occasionally. Further, a straightforward argument suggests that the problem should get 
worse as the lattice spacing is made smaller. In lattice units, in terms of the lattice inverse 
temperature the typical magnetic field strength is B ~ (3 L 1 , and since the violation 
of the Gauss law depends on the nonabelian nature of the theory the size of the induced 
violation must go as B x B ~ (3j}. To keep 11 fixed in physical units as we change the 
lattice spacing, \i should go as /3£ , so the induced unphysical E field goes as fxB x B cx 
and the spurious energy introduced per unit lattice 4- volume goes as /3^ 4 . The rate of Ncs 
violation per lattice 4-volume, for nf3 L fixed, due to genuine infrared topology changing 
processes should go as /^ 4 , and so the energy shift from these processes goes as (3j^ . If this 
naive dimensional argument is right then, if the Gauss constraint is unquenched or poorly 
quenched, the contribution to T M due to excitation of the unphysical modes should grow 
linearly with inverse lattice spacing, and no fine lattice spacing limit for T M will be found. 

The problem can be understood and cured by considering the implementation of the 
Hamiltonian equations of motion, Eq. (|28|), more carefully. When we modify the electric 
fields we should not modify all fields, but only the linear combinations orthogonal to those 
constrained to be zero. That is, the electric fields can be partitioned into two orthogonal 
subsets, the constraints, which I call the E c , and the fields orthogonal to the constraints, 
which I call the E* . At fixed U the E* are in correspondence with the momenta of the 
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cannonical basis of the observable subspace. Only these E* are dynamical, the E c should be 
held zero. Since U is fixed during the update of E in the leapfrog algorithm, this means we 
should modify the update rule in Eq. (|28|) to 

AE C AE C , 

— — = -—^ = 0+0 = 

AE* A E* 

-£f = -^(/x = 0) + /i{£*,iV C5 }. (31) 

This is exactly what we would conclude if we defined ANcs as E* ■ B/(2tt 2 ), which is 
equivalent to E ■ B/(2tt 2 ) on the constraint surface, since the E c are zero there. 
We can implement this update by finding an orthonormal basis for the E*, 

Ea = ^2 C aiEi , ^2 C ai c Pi = $af3 (32) 

i 

and updating them by changing them by 

K = K-^Y,** B i (33) 



or 

fiAt 



E i = E i~ 7^2 51 C d C OJ B i ' ( 34 ) 



27T 2 

a j 



which is basis independent and preserves the Gauss constraint by construction. Equivalently, 
given a complete orthonormal basis of Gauss constraints E c , 



E c — 2J d a iEi , daidpi — 5 a p c a iC a j + ^ dpidpj — 8ij , (35) 

% i « p 



the update is 



Ei = Ei ~ 2^ \B i -J2daiJ2 d ^ B j ] • l3Gi 



a j 

This technique as described is impractical. The problem is that (unlike in the abelian 
theory) the linear combinations E* are in general not well localized and change with each 
time step. Finding them is a problem in the diagonalization of sparse matricies, but even if 
it could be performed efficiently, the number of operations required to implement Eq. ( |3"1D 
grows with the square of the number of lattice sites, which makes it numerically impractical. 
We also cannot implement the algorithm by using Eq. (RB) because the natural, local basis 
for the Gauss constraints is not orthogonal, and again any orthogonal basis is not localized 
and changes with time. 

However, there is a (linear in the number of lattice points) algorithm for implementing 
this update scheme approximately, which can be made highly accurate. The idea is that any 
update of form 

Ei — > Ei — \ B i — ^2 daiK-aj , (37) 
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with d a i describing any basis for the Gauss constraints, and with the n a completely arbitrary, 
will correctly modify the E*. It is not even essential that the basis for the Gauss constraints 
be orthogonal, as long as it is complete, independent, and orthogonal to all E*. What 
uniquely determines the right choices for the K a (once we have chosen a basis for the Gauss 
constraints and hence the d a i) is that the Gauss constraints must be satisfied at the end. If 
we can find an algorithm which, by only making changes to the E fields orthogonal to the E*, 
incompletely but very accurately enforces the Gauss constraints at each step, by choosing 
almost the correct K a , then we will have almost the evolution of Eq. (|34D ; and by improving 
the accuracy to which the Gauss constraints are approached we can test whether the results 
depend on the residual failure. This is the technique I will adopt for the evolution of the Yang- 
Mills system with a chemical potential. It is still not clear whether the failure to correctly 
simulate the ultraviolet physics will influence the resulting T M (although I will of course 
check if my results have sensible small volume behavior), but the infrared evolution should 
be correct, and no spurious physics should arise from violations of the Gauss constraints. 

Restoring Gauss Constraints 

I will now briefly describe two algorithms for the removal of the Gauss constraint vio- 
lations. For the uninterested reader it is only necessary to know that algorithm 1 removes 
about half of the accumulated Gauss violation at each step, and algorithm 2 is about twice 
as costly in compute time, but leaves a residual Gauss constraint about one fifteenth as large. 

The basic idea of quenching the Gauss constraint is to change the E field on a link in the 
direction which will reduce the Gauss violation at both endpoints. The most basic method 
is 

Ei{x) -> Ei(x) + i{UiC{x + l)U} - C{x)) (38) 

where by C(x) I mean the Gauss violation at the point x, and I have of course parallel 
transported the forward C to the point where the group indicies of E{ reside. This algorithm 
is equivalent to the relaxation algorithm with Hamiltonian H = J2x C(x) ■ C (x) (Lie algebra 
dot product) used in |J, and it does not alter the E*. Algorithm 1 is to apply this relaxation 
once at each timestep, using a value of 7 chosen to make it efficient. To evaluate its efficiency 
it is useful to use Fourier analysis, which should be approximately valid in the ultraviolet 
(where most of the Gauss violation occurs). When C(x) = C(k) sm(k ■ x), each update takes 

C(ife)-> |l-7^2(l -cos ki)^jC(k). (39) 

The sum is recognizeable as u> 2 (k) for the lattice dispersion relation, and its maximum value 
is 12. The sign of C(k) oscillates with each application of the algorithm if 7 > l/u) 2 (k), and 
if 7 > 1/6 then the most ultraviolet mode becomes unstable. For 7 close to 1/6 the damping 
becomes inefficient for the most ultraviolet modes. However, the damping of infrared modes 
is in general quite inefficient, and it is good to have 7 as large as practical. Using 7 = 0.1 
or 0.12, I find the algorithm is quite efficient. If the algorithm is applied at each step, then 
the total Gauss violation, measured by \Jj2 x C( x ) ■ C(x), after a step is about as large as 
the Gauss violation generated by one step, starting from no violation. For an 18 3 grid, a 
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timestep of 0.06 lattice spacings, a chemical potential fi = 0.2 in lattice units, and a lattice 
inverse temperature (3^ = 6, I find (C(x) rms ) ~ 1.7 x 10~ 4 . 

To determine how much of the accumulated Ncs i n the realtime evolution of the lattice 
system was due to the Gauss constraint violation, I performed a linear regression fit of 
the energy of a lattice system to Ncs when the chemical potential term was turned on 
with /i = 0.2 in lattice units, at (3l = 5 (// = T). The slope should equal /i if no Gauss 
constraint violations are generated; it was low by 7.5%, indicating that a nontrivial minority 
of the measured change in N C s was due to the violation and subsequent quenching of Gauss 
constraints. When the algorithm was only applied every 16 steps, the slope was low by 37%. 
Clearly, then, quenching the Gauss constraints effectively makes a great deal of difference, 
but algorithm 1 is insufficiently thorough. 

Algorithm 2 is a modification and improvement of algorithm 1. Because the time stepsize 
is by necessity much shorter than any frequency in the problem, the Gauss constraint vio- 
lation generated at each time step will approximately equal that of the previous time step; 
a good first correction is to repeat some multiple m < 1 of the total modification made in 
the previous time step. Then one application of the relaxation algorithm described above 
is applied. This combination leaves an rms Gauss violation a factor of 4 smaller than algo- 
rithm 1. A further improvement is accomplished by applying the relaxation step twice each 
time step, with different stepsize constants 71 and 72. All that is necessary for algorithm 
stabiltity is that 1 > (1 — 7iu; 2 )(l — 72CU 2 ) > —1/3 (the lower bound would be -1 for m = 0, 
but changes when m^0. The value I use is for m = 1.). By making the 7 well separated, 
with 72 > 1/6, I make the algorithm efficient over a wide frequency range. I find that the 
parameter choice m = 0.9, 71 = 5/48, and 72 = 5/24 gives excellent performance, with the 
Gauss constraint violation after application about 1/15 as large as with algorithm 1. Using 
this combination, the linear regression fit of energy versus Ncs gives /1 to better than 1% for 
every lattice spacing used in this paper. I conclude that evolving the system using algorithm 
2 is for all practical purposes equivalent to implementing Eq. ([34]). Incidentally I can also 
conclude that I implemented the chemical potential with the right numerical coefficient. 

If it should prove necessary, it is straightforward to further improve algorithm 2; one 
repeats the relaxation step more than twice, with a different stepsize constant each time, 
chosen so that the function U(l — / -f i u) 2 ) will be as near zero in as wide a range as possible, and 
never outside the stability bound. For instance, I found that 3 appliacations with 7 = 7/72, 
7/48, and 7/24 leaves a residual violation about 2.5 times smaller than algorithm 2. Also, 
reducing the stepsize improves algorithm performance; algorithm 1 improves linearly with 
At and algorithm 2 (keeping (1 — m)/At fixed) improves quadratically. 

4 Some Numerical Results 

I implemented the Hamiltonian system with chemical potential described above, and the 
thermalization algorithm of Krasnitz |Tl| . The time stepsize for the Hamiltonian evolution 
was chosen as follows: in the ultraviolet (where most of the energy of the system resides) the 
Fourier modes behave as weakly coupled harmonic oscillators with frequency given by the 
lattice dispersion relation uj 2 = 2J2{^ — cosfcj). Defining the energy of the system at time t 
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to be 

Energy (t) = £ 1 - £lH7 D (f) + £ ^± 2 j + £ ^ aJ , (40) 

the leapfrog algorithm should keep the central value of energy stable for a harmonic oscillator, 
but the amplitude of oscillations in the energy of a harmonic oscillator should be u 2 At 2 /A 
times the energy. To keep these fluctuations < 1% for the most ultraviolet mode (uj 2 = 12) 
it is necessary to choose At < 0.06. I choose At = 0.06 in all simulations discussed in this 
paper. Integrating over all the modes, the perturbative estimate of the variance in the energy 
of an N 3 lattice at lattice inverse temperature /3l, in the limit of large N, is 



for my value of At. The energy of the system without chemical potential was in fact ab- 
solutely stable to overall drift, with variance within 10% of the analytic estimate, for runs 
on large and small lattices. I also checked that, for initial conditions with only one excited 
mode, the dispersion relations were correct. 

To test the thermalization algorithm I measured the overall energy and compared it to 
QN 3 /f3i, which is the free field estimate (there are 6 degrees of freedom per site); it agreed to 
a few percent, better at large (3l and worse at low (3l- I also compared the values of Wilson 



loops to the results of Krasnitz |TT[ , and found good agreement. 

A more substantive check was to measure the Ncs diffusion rate in the absence of a 
chemical potential for a thermalized system, comparing to the results of Ambjorn and Kras- 
nitz. Because my computer resources were limited I could not check their results at the 
large values of (3l they used; instead I measured Yd at a range of smaller (3l, to find at 
what value lattice coarseness effects arise. I was also unable to demonstrate in a con- 
vincing way that the ultraviolet, white noise fluctuations satisfy their analytic estimate of 
SN CS = 0.00684A^ 3 /(7r/?i) 2 ; instead I assumed this behavior and subtracted it off from the 
value of (Ncs(t) — N C s(fy) 2 ] however I compare two values of t to check the validity of this 
procedure. Expressing T d = /^/(/Jl 7 '") 4 , (this value of k corresponds to the value defined 
in the introduction), I present my values for in Table 1. All error bars are one a and 
statistical. 

It is clear from the table that values of (3l < 5 are contaminated with some finite lattice 
spacing effect, but the results above this scale are consistent with those of Ambjorn and 
Krasnitz. 

Next I implemented the chemical potential on the lattice, as described in the previous 
section. As a first check to see whether the technique will be plagued with ultraviolet diver- 
gences, I thermalized the system so that all U and E would fall in an abelian subspace of 
SU(2) and its Lie algebra. This is easy to do with the thermalization algorithm of Krasnitz, 
because initial conditions which only contain excitation in one Lie algebra direction never 
have the other Lie algebra directions excited by the algorithm. The evolution is then equiv- 
alent to compact U(l) gauge theory, in which Ncs should oscillate about but never drift. 
Indeed, when I applied a chemical potential, the value of N C s did not drift, but fluctuated 
in a narrow range about 0; so if there are ultraviolet problems in the chemical potential 
method, they only arise out of nonabelian interactions. 
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As a next test of the reliability of the algorithm I measured (Ncs) with // = 1/Pl, which 
should give half T d V. I will express my results as k(hPl) = (Ncs)(Pl^) A /{^PlN 3 ); we 
should expect k([i[3l) = Kd/2 for small \iPl- I find that, for (3l = 5, «(1.0) = .46 ± .03; 
for /5 L = 6, k(1.2) = .50 ± .04; and for /5 L = 8, «(1.0) = .60 ± .11. All three results are 
on lattices of length N = 3/3l, well larger than the size Ambjorn and Krasnitz found would 
remove finite scaling effects. These results are in very good agreement with the measured 
values of at the same values of Pl, which is suggestive that both techniques are working 
correctly. 

However, since the implementation of the chemical potential, while probably correct for 
infrared excitations, is almost certainly wrong for ultraviolet excitations at the lattice scale, 
it is a good idea to test the finite volume behavior of the results. We expect that, when 
the volume of the lattice is made small, there will not be enough room for energetically 
unsuppressed topology changing events, which must be spatially extensive, to exist, and the 
rate of topology change should fall. If our results arise from true infrared behavior then we 
should see k(/j,Pl) fall; but if our results arise from spurious ultraviolet effects, then the rate 
should simply scale with the 4- volume of the simulation and k(/j,/3l) should be unchanged. 

I measured k(ixPl) for \ifi L = 1.2 and (3l — 6 for a number of lattice sizes. The results 
are presented in Table 2 and Figure |j. Each datapoint represents a total lattice 4-volume 
of about 2 x 10 7 a 4 (a the lattice length). It is clear that the rate is falling for small lattice 
size, and the curve is about the same as that found by Ambjorn and Krasnitz for Kd- Hence 
we can conclude that the results are in fact due to infrared physics, which the model should 
treat properly. 

With the algorithm thus tested, I can now investigate a new question; what is the behavior 
of the system under large chemical potentials \x 3> 1/Pl? This directly tests whether there is 
a free energy barrier to the flow of Ncs, as 1 discussed in section 2. Again, I use Pl = 6 in a 
compromise between compute time (which goes as f3\) and lattice refinement. I present some 
results in Table 3 and plot them in Figure [5| For the larger values of \x it was necessary to 
correct for the finite heating of the system, and to use several short evolutions of independent 
thermal initial conditions to minimize the amount of heating. The data show that (Ncs) 
rises almost linearly with fi up to /i/^ > 6, in stark contrast to the behavior when there 
is a free energy barrier. The errors quoted are all statistical, but there may be common 
systematic errors at the level of 3% from thermalization and finite stepsize. I also measured 
(N cs ) for p L = 10, nP L = 10, and iV = 30 and found k(10) = .607 ± .022, which shows 
that the approximate linearity is not a small Pi artefact, but that the slight departure from 
linearity (the excess of k(10) over k(1)), is also probably real. 

5 Modified Lattice Hamiltonians 

There is one potential problem with the classical lattice technique which might mean that the 
results for T d and do not correspond to the correct continuum, quantum behavior. While 
the thermodynamic properties of the lattice system may be the same as the continuum 
theory in the infrared (in the very good approximation of dimensional reduction for the 
continuum theory), the dynamics may be different; in particular, while we know that the 
"hard thermal loop" corrections to the static propagator, namely the Debye screening, are 
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of the same functional form as in the quantum system, this does not extend to the nonzero 
frequency case. Bodekeret. al. have recently shown that, because the most ultraviolet modes 
move under lattice rather than ultrarelativistic dispersion relations and have a rotationally 
noninvariant ultraviolet cutoff, they induce corrections in the propagators of the intermediate 
scale which are not the same as the "hard thermal loop" effects of the quantum theory, and 
are in fact not even rotationally invariant |2l|0. This sort of cutoff dependence is typical of 
linearly divergent quantities. While the lattice coarseness independence of demonstrates 
that the overall magnitude of the "hard thermal loop" effects does not matter (at least when 
the lattice is fine enough that the Debye screening scale and the nonperturbative scale are well 
separated), it is possible that the detailed k/ko dependence of these corrections matters in 
the determination of rV This possibility is motivated by Braaten and Pisarski's calculation 
of the gluon damping rate, in which the magnitude of the plasmon mass (proportional to 
the strength of hard thermal loop effects) cancels out, but the functional form of the hard 
thermal loops determines the calculation [p^j . 

A related but less urgent matter is to understand why T d falls as the lattice becomes 
quite coarse. There are two possible explanations for this behavior. One is that the lattice 
is incorrectly evolving the infrared equations of motion due to artefacts, namely nonrenor- 
malizeable operators induced by the crude choice of lattice action. As the lattice becomes 
finer the effect of these artefacts on the behavior of the nonperturbative infrared modes falls 
as (5J 2 . Another possibility is that Td is only constant in the limit that the nonperturbative 

1 /2 

length scale [0(/3l) in lattice units) and the Debye screening length scale (0(/? L ) in lattice 
units) are well separated. 

One can study at least the first of these issues by the following technique. I design a 
lattice Hamiltonian which produces different (wrong) ultraviolet dispersion relations and 
will induce different (wrong) hard thermal loops. Then I measure for different lattice 
coarsenesses and see if there is a large lattice size limit and whether it is the same as the 
limit for the normal cubic lattice action. If it is, then the limit did not depend on which 
(wrong) hard thermal loops were induced, and would presumably not change if the right 
hard thermal loops could be induced. 

I examined two modifications of the lattice action. In the first, I make the spatial part of 
the Hamiltonian 5/3X)o(l — l/2TrZ7n) — 1/12 X^(l — l/2Trt/| = |), where the second sum is 
over all 1 x 2 rectangular plaquettes (both orientations). This choice gives steeper dispersion 
relations in the ultraviolet. For modes with d\Ei = d^E 2 = d^E 3 = the expansion of 
u 2 in k has no quartic term; but for other modes it has positive k\k\ type terms, and the 
action is not an "improved" action in the sense of eliminating dimension 6 operators from 
the Hamiltonianf]. 

I analyzed this model for several lattice coarsenesses (values of /3l), generating the initial 
conditions with the thermalization algorithm discussed in Appendix A. Because the Hamil- 
tonian is more complicated, the numerics were more time consuming; to get good statistics 



6 In fact Bodeker et. al. show this explicitly only for the abelian Higgs theory but there is no reason to 
doubt that the same thing happens in Yang-Mills theory as well. 

7 The easiest way to see this is to consider the thermodynamics of the model; introduce Lagrange mul- 
tipliers Aq(x) for the constraint at site x and perform the Gaussian integration over the fields E. This 
generates a kinetic term for A which is precisely the minimal (unimproved) implementation of (DiA ) 2 /2. 
An improved action should produce an improved kinetic term for A . 
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I was forced to assume that (Ncs) rises linearly with /i up to (3 fx ~ 5, rather than verifying 
it explicitly as I did for the standard Hamiltonian. Some results are presented in Table 4. 
The same scaling behavior is reached, although the lattice must be finer before it is attained. 
It is not clear whether the approach to scaling is slower because the Debye mass is smaller 
at the same (3l (it depends on the integral of u~ 2 , and u is larger in the ultraviolet in this 
model), or because of the different dimension 6 artefacts. 

An alternate and stronger check is to deliberately make cutoff scale lattice behavior even 
worse, so the hard thermal loops will look even less like their correct values. If the small 
lattice size limit is unaffected then there can be little doubt that the functional form of the 
hard thermal loops really is not important. 

To do this I make one lattice direction longer than the other two. That is, defining the 
electric field E by U(t + At) = exp(iAtr ■ E)U(t), I take as my dot product 



A-B='£l(A 1 B 1 + A 2 B 2 ) + r 1 A 3 B 3 , (42) 



sites 



as my Hamiltonian 



H = ^ + E *(1 " ^TrC/ Dl2 ) + r\l - l -TvUn 13 + 1 - \^Un 23 ) , (43) 



as my magnetic fields 



85? = n^TrHr")^, 

8D Z 

m = l-^\^{-ir a )U U31 , 

8B« = l^Tr(-iT a )U Dl2 , (44) 



as my definition of DjF 3l 



(DjFj^ = l-^\^{^ a )U nri +Y,\^(-^)Uu l2 

□13 Z Ol2 Z 

(DjF j2 ) a = l-^lM-^ a )U D23 +J2l-TT(-tT a )U Dl2 

□23 A Oi2 Z 

{D^) a = ^ilVHrJC/a, (45) 

(where the orientations should be taken as appropriate), and as my AN C s 

E-BAt , x 

AN CS = - 2 . (46) 

271 

The update rules for U and E are, in terms of B and DjFji defined above, unchanged. The 
thermalization algorithm and the Gauss law restoration are discussed in Appendix A. The 
thermalization is applied so that the total energy will be 6lN 3 /(3 (6 the number of degrees 
of freedom). 
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To determine how we should expect this modification to change the observed rate of Ncs 
violation, it is simplest to successively stretch each lattice direction by the same factor of Z. 
The resulting model is the same as the original cubic lattice model, except for factors of Z. 
The energy in the magnetic field goes as Z -1 J2n 1 — l/2TrZ7 n , but its magnitude grows as Z 3 ; 
so the magnitude of 1 — l/2TrZ7n has grown as Z 4 . Therefore the new model corresponds to 
the old model with @ L divided by Z 4 . Also, because the expectation value of E only grows as I 
and not Z 2 , we find that the time step has changed by a factor of Z _1 . The chemical potential, 
in lattice units, is unchanged. Therefore, the expected value of (Ncs) should change as the 
4*3 — 1 = 11 power of Z, and stretching only one direction by Z, using the definition of the 
system discussed above, should change the rate by Z 11//3 . To compensate, I shift my value of 
Pl by Z _4//3 and my value of At by Z 1//3 ; I should then expect the rate to be Z independent, 
as long as /3l is large enough to eliminate finite spacing artefacts and iV is large enough to 
eliminate finite volume artefacts. Because the grid is rectangular, these are more stringent 
conditions than in the cubic lattice case. 

I have tested the behavior of this streched lattice for several values of Z, choosing /3l 
(including the Z~ 4 ^ 3 correction) to be around 8. Some results are presented in Table [8[ 
The data suggests a very weak dependence on Z, with the rate rising slightly for quite high 
values. This suggests some weak lattice shape dependence, presumably arising from a weak 
dependence on the hard thermal loops, but the effect is small and it is reasonable to believe 
based on this evidence that the value for k determined from the cubic lattice should be quite 
close to the value we would find if we could include the right hard thermal loops. 



6 SU(3) 

The topological structure which exists in SU(2) (weak isospin) gauge theory also exists 
for SU(3) (color); there is a Chern-Simons number, and at high temperatures, thermal 



excitations should be the dominant means of topology change [p3j] . The rate of these changes, 
Tstrong, is another important quantity to baryogenesis, since it may control the rate at which 
fermions change chirality in the plasma pij. Parametrically, the Ncs diffusion rate can be 
written as r dstrong = K d)Strong (a; s T) 4 , with the constant n dstrong another unknown. 

The rate can be measured by a straightforward extension of the chemical potential 
method used for SU(2). Group and Lie algebra elements should be chosen from SU(3) 
rather than SU(2), and 1 — l/2Tr should be replaced by 3/2 — l/2ReTr; with these sub- 
stitutions everything carries over. The lattice inverse temperature is now related to the 
physical temperature as /3l = 4/(g 2 aT), which will generate the correct power of g s ; the 
constant /^strong = 2K MiStrong can be gotten from simulations in the same way as k^. The 
thermalization algorithm, presented in Appendix A. also carries over to SU(3) without seri- 
ous modification. 

In practice, because the group SU(3) is larger than SU(2), and because the anticom- 
mutators of the Lie algebra generators A a are not multiples of the identity element, it is 
much more numerically costly to simulate the group SU(3); I find that the time to update 
a lattice of the same size is approximately 8 times larger for SU(3) than for SU(2). Also, I 
find that the inverse lattice spacing (3i must be larger before the results become (3i indepen- 
dent. Because of these complications, I was unable to verify that Ncs rises linearly with \i. 
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Instead, I have assumed that this behavior carries over from the SU(2) case. I can then use 
quite a large value for /j,/3l, which is necessary to get good statistics out of relatively short 
runs. (The length of a run necessary to get fixed statistical error scales as the square of fi.) 
Unfortunately, this introduces a possible systematic error; it can be eliminated only with a 
great deal more computer time. 

I have also been unable to check explicitly that the rate displays similar lattice volume 
behavior to the SU(2) case; instead I have assumed that the rate becomes volume indepen- 
dent on lattices with N > 2@l. Since the magnetic screening length of SU(3) is presumably 
shorter than that for SU(2), this assumption is almost certainly justified. 

I present some numerical results for /^strong in Table [|. Because of the limitations dis- 
cussed above, these should be considered preliminary; the systematic error bars are probably 
at least 10%. A preliminary estimate for the ratio /^strong / K n is 7 ± 1, but it is not clear 
from the data whether the finest lattice used was large enough to reach the fine lattice limit. 

Note also that the energy is systematically higher than the leading order perturbative 
estimate (the number of degrees of freedom divided by (3l). This effect is real. Since the 
electric fields appear quadratically in the Hamiltonian, I can compute that the energy in E 
fields should be 8N 3 //3l, and this is satisfied numerically. The same energy excess in the 
magnetic fields was observed in SU(2); for SU(3) it is more than twice as large, presumably 
because there are more ways for the fields to interact. 

7 Conclusion 

It is possible to include a chemical potential for Ncs in lattice simulations of the classical 
Yang-Mills field equations, and the rate of Ncs motion is related to that of Ncs diffusion 
in the manner demanded by the detailed balance argument. Furthermore, the results show 
that there is no free energy barrier to topology change in the symmetric phase, because the 
response of Ncs grows linearly with chemical potential over a much wider range than is 
possible if the transitions occur by hopping over a barrier. I have also demonstrated that 
the classical, lattice rate is weakly, if at all, dependent on the details of the hard thermal 
loops induced by the high momentum modes. 

To conclude, I will mention a possible interpretation of the results, and some directions 
for future work. 

The two striking things about the rate of Ncs change under a chemical potential are 
that, in terms of the volume necessary to remove finite size effects, the rate seems very low, 
and that the rate remains linear up to surprisingly large values of chemical potential \i. 
These two things may be related. Consider a square Wilson loop of length x on a side in the 
continuum limit of the classical theory. When x « (3 (by which I mean a/?£ as we take (3^ 
to infinity and a to zero) the trace of the loop is 2(1 — Cx/(2j3)), where 

n 2 - 1 r 1 - cos(fci) 1 - cosffe) r, n ; , 9 
C = — — / dk x dk 2 rj^V^i + fc 2 (47) 

for SU(n) gauge theory. This expression has an ultraviolet logarithmic divergence, but what 
I am interested in is the contributions from k ~ ir/a, which are 0(1), with no powers of tt. 
The length scale where the "correction" term which makes the trace ^ 2 becomes order 1, 
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which is the scale where nonperturbative physics sets in, is x ~ (3l, with no powers of it. This 
is the length scale where excitations no longer remain coherent long enough for energy to 
oscillate between E and B fields, and it is reasonable to expect, then, that it is the scale where 
Ncs will lose its oscillatory behavior and behave in a Brownian fashion. The characteristic 
4-volume of Brownian motion is thus 0(/3 4 ), with no powers of n, which corresponds to the 
minimum volume for finite volume effects to disappear, and the magnitude of E-B, integrated 
over one such 4-volume, is 0(1), again with no powers of 7r. To get SNqs = E ■ B/(27r 2 ) 
to be 0(1) then requires 0(7r 4 ) characteristic Brownian motion 4- volumes, and we should 
anticipate T ~ (7T/3) -4 , not /3~ 4 , or ~ 1, not ~ 7r 4 . 

Turning to the chemical potential, one would expect that \x would start seriously changing 
the dynamics of the infrared when [iE-B /(27r 2 ), integrated over a volume of order /3 4 , is order 
that is, when nf3 ~ ir 2 ; and that for values smaller than this, the chemical potential is a 
small correction to the dynamics and the response should be linear. Again, this corresponds 
well to the actual behavior of the system. 

What possible future projects can the chemical potential method be applied to? One 
question I could not answer here was whether the Ncs diffusion rate falls off on coarse lattices 
because the Debye screening length comes on order the scale of nonperturbative infrared 
physics, or because of the effects of nonrenormalizeable operators arising from the crudeness 
of the lattice action. One could resolve this question by finding a lattice Hamiltonian system 
and an update algorithm which do not generate dimension 6 nonrenormalizeable operators 
in the infrared physics. If dimension 6 nonrenormalizeable operators are causing problems 
on coarse lattices then such a program could also allow better study of the SU(3) case 
without requiring the very fine lattices which make such study numerically expensive. Also, 
knowing whether or not Debye screening is important would help to clarify the role which 
hard thermal loops play in the calculation. 

Also, extending the "chemical potential method" to the Yang-Mills Higgs theory allows 
one to study the out of equilibrium decay of the infrared gauge fields as they are swept 
up a bubble wall during the electroweak phase transition. It is frequently assumed in the 
literature that any left handed baryon number asymmetry generated by the motion of the 
bubble wall must find its way to the symmetric phase to cause baryon number generation, 
because N C s violating processes are exponentially suppressed in almost all of the wall and in 



the broken phase[25|; and a contrary assumption, that a chemical potential for Ncs present 



on the wall can act on the infrared gauge field configurations as they decayf26|, has never 



been convincingly tested. The method proposed here could be used to do so, for instance by 
integrating out the fermions analytically to produce effective interactions in the two doublet 
model, and numerically evolving the system through the electroweak phase transition. 
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8 Appendix A 



In this appendix I present a simple, efficient thermalization algorithm for the classical lattice 
Yang-Mills or Yang-Mills Higgs system, or for any system in which the phase space is a vector 
bundle over the space of coordinates, with a Hamiltonian which is a pure coordinate term 
plus an inner product on the vector space, with the inner product allowed to be coordinate 
dependent. Constraints on the momenta are permitted provided they are linear, ie in each 
vector space they restrict to a vector subspace. (In the case of Yang-Mills theory the phase 
space is the tangent bundle over the space of U, and the constraints restrict the momentum 
vector to have a zero projection in each direction which corresponds to a gauge transformation 
on the space of U.) The algorithm is faster than that of Krasnitz, and much easier to apply to 
new systems, although the above requirements make it less general. Almost all the numerical 
equipment it uses must be developed anyway for the evolution of the system with a chemical 
potential. 

The idea is the following; we begin in the vacuum. Then we fully thermalize the E fields, 
leaving the link matricies U untouched. That is, we draw E with correct Boltzmann weight 
from the vector space of allowed E at that value of U. We then evolve the fields according to 
the classical Hamiltonian equations of motion, allowing the thermalization to mix between 
E and U fields. Then we throw out the values of the E fields and fully rethermalize them; 
we repeat the thermalization and evolution over and over until the system is completely 
thermalized. 

To see that this is a correct thermalization algorithm, consider the evolution under this 
algorithm, averaging over realizations, of an ensamble of systems described by some proba- 
bility distribution. If we really draw E with correct thermal weight from the fixed U fiber 
(the set of E at fixed U), then the correct thermal probability distribution will be preserved 
by the algorithm. Further, any probability distribution preserved by the algorithm must 
have the property that the probability distribution is a product of a function on the coor- 
dinate manifold and a function on each fiber, with the function on each fiber the thermal 
distribution for the E fields at that value of U ; and this property must be preserved by 
the Hamiltonian evolution. Only the correct thermal distribution satisfies these conditions. 
Hence the algorithm should be correct. The only problem is finding a way to draw E from 
the thermal distribution on a fiber. 

Now the thermal distribution function on the fiber of a particular U (that is, on the 
space of E at fixed U) is Gaussian. All we have to do to thermalize the E fields in this 
vector space is draw a value from the Gaussian distribution. In the absence of constraints, 
this would be easy. The problem is that, in the obvious natural basis for the E fields, the 
constraints are nontrivial linear combinations of E fields, and to find an orthonormal basis 
of the directions in the dynamical subspace (the subspace permitted under the constraints) 
we must diagonalize a large matrix, as discussed in section 3. However, this does not turn 
out to be a fundamental problem. 

The key is to notice that the obvious basis of E fields is orthonormal (up to a factor of 2) 
under the inner product defined by the Hamiltonian. If we choose each E field independently 
from the Gaussian distribution, then that will also choose the fields independently from the 
Gaussian distribution in any other orthonormal basis for the E fields, a special property of 
the Gaussian distribution. This includes any basis which decomposes into a subset {E*} 
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orthogonal to the constraints and a subset {E c } along constraint directions. This means 
that we will have chosen the E* correctly. All that remains is to set the E c to zero without 
changing the E*, that is, to orthogonally project to the constraint subspace. But I have 
already presented a dissipative algorithm for doing so; by applying the dissipation step 
repeatedly, we can quite easily bring the magnitudes of the Gauss constraint violations to 
the level of machine roundoff accuracy. The electric fields are then correctly thermalized. 

I should comment that the demands on the dissipative algorithm here are different from 
the demands during the evolution of the fields with chemical potential; so the optimal values 
of the algorithm coefficients are quite different. I find that, for algorithm 2, the value of m 
should be brought down to m < 0.5. 

I implemented this algorithm for the Yang-Mills system and compared it to the algorithm 



of Krasnitz||ll||; I used the stepsize At = 0.03 in lattice units for my algorithm and At = 7^ = 
0.03 for the algorithm of Krasnitz. I used 10 E field randomizations with evolutions of length 
/3l between them for my algorithm. On lattices of iV = 15 and N = 18 I systematically 
found my algorithm gave total energies about 1% less than the algorithm of Krasnitz; this is 
probably because of the rather large time step and dissipation constant I used to implement 
his algorithm. (In his paper he advocates a more conservative stepsize.) I also compared 
the values of Wilson loops and found them to be within the same 1% margin. That their 
values are preserved by the Hamiltonian evolution of the thermalized system, is more or less 
automatic in my algorithm; but I tested it anyway and found it was true. 

It is completely straightforward to extend this algorithm to the group SU(3); my results 
for SU(3) use my thermalization algorithm. It is also reasonably straightforward to apply 
it to SU(2) Higgs theory; but here it must be remembered that the constraint is a linear 
combination both of E fields and of n fields (the conjugate momentum of the Higgs doublet); 
the dissipation algorithm used in the thermalization should change both fields in proportion 
to their contribution to the Gauss constriant (holding U and <fi fixed throughout the enforce- 
ment of the constraint). The same algorithm must be used in the evolution of the system 
under a chemical potential. 

It is never necessary to dissipatively enforce the Gauss constraint during the evolution step 
of the thermalization algorithm, because no chemical potential is applied and the equations 
of motion preserve the constraint identically, as I have of course verified numerically, for 
cubic and distorted lattices and for SU(2) and SU(3). 

When applying the algorithm to the distorted lattice, it must be remembered to define 
orthogonality in terms of the quadratic part of the Hamiltonian. In other words, since 
Energy=Z(£' 2 + E\ + (Z -1 ^) 2 ), when we find the violation of the Gauss constraint C x = 
Ei + E 2 + l~ l (l~ l Ez) — (incoming lines), we should apply a correction of —^C x to both Ei, 
and E 2 , and of — 7(/~ 1 C :c ) to (Z -1 ^); so although l~ 2 E 3 enters the Gauss constraint, E 3 
is modified by — 7C X , with no power of I. The same dissipation algorithm should be used 
during the evolution of the fields with chemical potential as in the dissipative part of the 
thermalization, because the inner product defined by the dot product is the same as that 
defined by the Hamiltonian. 

As a check that the algorithm correctly chooses E fields for the elongated lattice system, 
I thermalized the system with I = 1.5, ending by randomizing the E fields. I then allowed 
the system to evolve under the Hamiltonian equations of motion discussed in section 5, and 
checked to see if the distribution of energy between E ly E 2 and E 3 changed; it did not, to 
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within statistical errors of less than 1%, showing that the thermalization correctly partitioned 
energy between these fields. (If I shift E 3 by —l~ 2 C x in the dissipation algorithm, which 
would be appropriate if I used the obvious rather than the correct metric on the E fields, 
then the energy distribution does shift.) The relation between Wilson loops in the normal 
and elongated directions also was not affected by randomizing the E fields and allowing the 
system to evolve. 

I should comment that the thermalization problem cannot be solved by the following, 
very simple (but wrong), algorithm. Since the plaquette part of the Hamiltonian does not 
depend on the E fields, and the E field part is quadratic but with constraints, one uses a 
conventional lattice gauge theory algorithm to thermalize the plaquette part and then the 
above, simple algorithm to thermalize the E fields. The problem is that the constraints 
depend on the U fields; if we integrate over the (Gaussian) E fields we get a residual which 
is U field dependent and shifts the relative weights of different U field configurations; in 
fact it is exactly the A fields of the dimensionally reduced theory. What we could do 
is use a standard lattice thermalization algorithm to thermalize the 3 dimensional system 
with plaquette action and A Q fields with zero (bare) mass, use the resulting U fields, throw 
away the A fields, and choose the E fields as discussed above. This also ilustrates the 
relation between the thermalization algorithm presented here and the "molecular dynamics" 
Monte-Carlo algorithm of lattice gauge theory. 
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Figure 1: The 8 plaquettes contribute to the magnetic field along the vertical link. 



Pl 


3 


4 


5 


6 


7 


8 


Kd,t=/3 L 7T 


.48 ± .04 


.71 ± .06 


.93 ± .09 


1.16 ± .10 


1.07 ± .10 


1.14 ± .12 


Kd,t=2f3 L n 


.49 ± .06 


.67 ± .08 


.93 ± .13 


1.08 ± .13 


1.04 ± .14 


1.30 ±.20 



Table 1: Values of K d , the dimensionless T d , for various inverse lattice spacings. All data- 
points are for lattices of size N = 3(3l to prevent finite volume effects. Finite spacing effects 
become important around (5l — 5. 
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Link values at one time step Link values at P revious timeste P 

Figure 2: The configuration of non-identity links which will make the discrete spacetime 
integral of FF nonzero; Ul — 1 + ier-i, U2 — 1 + ier 2 , and U3 = 1 + ier 3 . Going around 
the plaquette containing Ul and U2 picks up an 0(e 2 ) contribution in the r 3 Lie algebra 
direction, which gives a nonzero contribution when traced against the plaquette which goes 
backwards in time, even though all links at large distances are the identity. 



N 


N/Pl 


k(1.2) 


4 


.67 


.15 ± .04 


6 


1.0 


.19 ± .05 


8 


1.33 


.41 ± .07 


12 


2.0 


.49 ± .06 


18 


3.0 


.50 ± .04 


24 


4.0 


.53 ± .06 



Table 2: Dependence of n(n(3 L ) on lattice volume for (3 L = 6. The rate of topology change is 
clearly suppressed on small volumes and saturates around iV = 2/3 L . Error bars are statistical 
and la. 
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Figure 3: Contributions to Z^i^ at a lattice point arise from the surface of the smallest box 
enclosing that point. Each link on the surface is traversed once with each orientation; in the 
abelian theory their contributions cancel, but in SU(2) theory, nonzero commutators arise. 
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N/beta L 

Figure 4: Lattice size dependence of the rate Ncs- The horizontal scale is N/Pl, the vertical 
scale is k(/i(3l)- If the rate arose out of ultraviolet physics then the curve would be flat. 



N 




(N cs )((3 L n) 4 /N 3 




18 


0.5 


.29 ± .08 


.58 ± .16 


18 


1.2 


.60 ± .05 


.50 ± .04 


18 


2 


1.15 ± .09 


.57 ± .05 


18 


3 


1.34 ± .09 


.45 ± .03 


18 


4 


2.01 ± .09 


.50 ± .02 


18 


5 


2.53 ± .09 


.51 ± .02 


18 


6 


2.92 ± .13 


.49 ± .02 


18 


7 


3.66 ± .13 


.52 ± .02 


18 


8 


4.19 ± .13 


.52 ± .02 


18 


9 


4.81 ± .14 


.54 ± .02 


18 


15 


8.93 ± .30 


.60 ± .02 


18 


16 


9.76 ± .40 


.61 ± .02 


18 


18 


10.0 ± .4 


.56 ± .02 


18 


20 


11.38 ± .46 


.57 ± .02 


18 


22 


13.1 ± .5 


.60 ± .02 


18 


24 


14.9 ± .6 


.62 ± .02 



Table 3: Rate of N C s change as a function of chemical potential. The rate is surprisingly 
linear, but turns up somewhat for /i(3 L > 10. Error bars are statistical and la. 




6 8 10 12 14 16 18 20 22 24 



chemical potential * bets^ 

Figure 5: Dependence of N C s on /i. Horizontal scale is fx(3 L , vertical scale is N C s(k(3 l ) 4 /N 3 . 
The line is a reasonable best fit to the points at lower \x. The dependence is startlingly linear 
up to ji ~ demonstrating that there is not a free energy barrier to change in N C s- 



Hamiltonian 


Pl 


N 




K 


new 


3.0 


9 


3.0 


.177 ± .011 


new 


4.0 


12 


3.0 


.281 ± .021 


new 


5.0 


15 


3.0 


.302 ± .032 


new 


6.0 


18 


3.0 


.357 ± .030 


new 


8.0 


24 


5.0 


.393 ± .034 


new 


10.0 


30 


10.0 


.585 ± .022 


old 


10.0 


30 


10.0 


.607 ± .022 



Table 4: Nqs violation rate for Hamiltonian including rectangular plaquettes. All volumes 
are large enough to eliminate finite volume effects. The rate begins to scale with lattice 
coarseness more slowly reaches the same limit, as shown in the last column. 
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A 


N 


I 






5.45 


24 


1.33 


4 


.56±.03 


6.70 


24 


1.14 


4 


.50±.04 


8.00 


24 


1.50 


6 


.73±.05 


8.00 


24 


1.50 


3 


.67±.05 


8.00 


24 


1.40 


6 


.64±.04 


8.00 


24 


1.30 


6 


.64±.04 


8.00 


24 


1.00 


4 


.57±.05 


8.00 


24 


0.80 


6 


.57±.04 


8.00 


24 


0.70 


6 


.54±.04 


8.00 


24 


0.60 


6 


.55±.04 


9.36 


24 


0.88 


4 


.65±.07 


10.77 


24 


0.80 


4 


.47±.08 


12.23 


24 


0.73 


4 


.60±.10 



Table 5: An assortment of data on an elongated lattice. There is no strong trend, although 
the very large I points turn up. 



Pl 


N 


Energy/(16iV 3 //3 L ) 




^strong (aO 


3 


9 


1.09 


3 


0.52±.05 


4 


12 


1.10 


4 


1.19±.06 


5 


15 


1.09 


5 


1.87±.09 


6 


18 


1.07 


6 


2.30±.ll 


7 


18 


1.05 


5 


3.01±.19 


10 


24 


1.035 


6 


3.19±.25 


12 


24 


1.032 


8 


3.86±.30 


14 


28 


1.025 


6 


4.26±.45 



Table 6: Rate of topology change in SU(3) lattice gauge theory at various lattice refinements. 
It is not clear that the fine lattice limit is reached, even at twice the refinement necessary in 
the SU(2) case. 
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